Alpha Decay in the Complex Energy Shell Model 
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Background: Alpha emission from a nucleus is a fundamental decay process in which the alpha particle formed 
inside the nucleus tunnels out through the potential barrier. 

Purpose: We describe alpha decay of ^^^Po and ^"''Te by means of the configuration interaction approach. 

Method: To compute the preformation factor and penetrability, we use the complex-energy shell model with a 
separable T=l interaction. The single-particle space is expanded in a Woods-Saxon basis that consists of bound 
and unbound resonant states. Special attention is paid to the treatment of the norm kernel appearing in the 
definition of the formation amplitude that guarantees the normalization of the channel function. 

Results: Without explicitly considering the alpha-cluster component in the wave function of the parent nucleus, 
we reproduce the experimental alpha-decay width of ^^^Po and predict an upper limit of T-1/2 = 5.5 x 10~^sec 
for the half-life of "*Te. 

Conclusions: The complex-energy shell model in a large valence configuration space is capable of providing a 
microscopic description of the alpha decay of heavy nuclei having two valence protons and two valence neutrons 
outside the doubly magic core. The inclusion of proton-neutron interaction between the valence nucleons is likely 
to shorten the predicted half-live of ^''*Te. 



PACS numbers: 23.60.-(-e,21.60.Cs,21.10.Tg,27.60.-|-j,27.80.- 



I. INTRODUCTION 

According to Gamow theory of alpha decay [l|, 0] , this 
fundamental radioactive decay can be considered as a 
two-step process d-Q. In the first step, an alpha cluster 
is formed inside the parent nucleus. The resulting alpha 
particle resides in a metastable state of an average po- 
tential of the daughter system. In the second step, the 
particle tunnels through the potential barrier. Each step 
requires different theoretical treatment. To compute the 
preformation factor that describes the alpha formation 
probability, one needs to evaluate the overlap integral in- 
volving wave functions of the parent and daughter nuclei, 
and that of the alpha particle. The estimate of the pen- 
etration probability requires a careful treatment of the 
resonance state. 

The commonly used formulation of the alpha-decay 
problem employs the R-matrix expression Q 



(1) 



for the absolute width. In this formalism, the first stage 
(formation of alpha particle with angular momentum L) 
is given by the reduced width 7^, while the second stage 
(decay) is expressed by means of the penetrability Pl- 
Alternatively, the absolute width can be obtained from 
the general reaction-theory expression [1-[IH 

Tl = Sl^I, (2) 

where Sl is the alpha-spectroscopic factor and T^^ is the 
single-particle (s.p.) decay width. 



Historically, expression ([T|) was derived in 1954 by 
Thomas Q using the time-independent R-matrix theory 
of nuclear reactions. In 1957, Mang [71 developed the 
alpha-decay formalism based on the time-dependent per- 
turbation theory. He made the connection with the shell 
model and succeeded in expressing the alpha-decay for- 
mation amplitude in a basis of s.p. states. As shown 
in Refs. [l3| formulations of Thomas and Mang are 
formally equivalent; there are, however, many differences 
when it comes to practical implementations. 

The reduced width calculated in a shell-model config- 
uration expressed in the harmonic oscillator (h.o.) basis 
is too small. This can be partly cured by means of con- 
figuration mixing involving extended shell-model spaces 
[ij, [Tsj as each admixed configuration contributes coher- 
ently to 7£. To improve asymptotic properties of s.p. 
wave functions, the particle continuum was taken into 
account by considering h.o. expansion [l7j or within 
a Woods-Saxon (WS) basis consisting of bound and out- 
going single-particle resonant (Gamow) states llSl. Il9ll . 
The configuration mixing calculations of Refs. [T^ [3 
in the valence space of ■^^^Po assumed the seniority-zero 
(pairing vibrational) wave functions obtained by consid- 
ering the monopole pairing interaction between like nu- 
cleons. However, all these improvements were not suffi- 
cient to reproduce the experimental alpha decay in ^^^Po. 
It is only after the valence proton-neutron interaction 
had been considered together with a generalized wave 
function expressed as a combination of cluster and shell 
model components [20l| that theoretical and experimental 
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widths could be reconciled [2lj. 

The R-matrix expression for the width ^ depends on 
the channel radius R. This radius should be chosen large 
enough so that the alpha-daughter interaction in the ex- 
ternal region is given by the Coulomb force alone [2^ . 
The infinite range of the Coulomb force implies, however, 
that the asymptotic behavior of the R-matrix expression 
is reached only at large values of R, at which the asymp- 
totic behavior of the shell-model s.p. basis (h.o. basis in 
most applications) used to calculate 7£(i?), does matter. 
Due to the mismatch between the internal part of the s.p. 
wave function (well described in the h.o. basis) and the 
asymptotic part (poorly or not described in the h.o. ba- 
sis), rather small changes in R may produce appreciable 
variations in penetrability. Physically, the reason for this 
sensitivity is the fact that the alpha cluster is formed in 
the surface region of the nucleus in which the coupling to 
the alpha continuum that impacts the radial behavior of 
the formation amplitude is important [16i] . Consequently, 
the absolute R-matrix width depends in general on the 
channel radius [2^ , and this is an obvious drawback 
of the method [24I . 

Our renewed interest in the alpha-decay problem is 
stimulated by the recent experimental data above the 
doubly-magic ^°°Sn [2^ that demonstrate the pres- 
ence of very fast alpha decays. Indeed, the observed en- 
hancement of the reduced widths of i05,i06rpg relative to 
2i3,2i2pQ jg two-to-three, thus confirming earlier expec- 
tations [26j of "superallowed" alpha decays in this region 
due to the large overlaps of valence s.p. shell model pro- 
ton and neutron wave functions. Our long-term goal is 
to estimate alpha preformation factors in nuclei above 
^°^Pb and ^™Sn by using large valence s.p. spaces, in- 
cluding positive-e nerg y Gamow states of a finite-depth 
WS potential [H, Il9j In this study, we focus on ^^^Po 
and ^°^Te nuclei having two valence protons and two va- 
lence neutrons outside doubly-magic cores. 

Our paper is organized as follows. Section |TT] briefly 
describes the alpha-decay formalism used in this work, 
with special emphasis on approximations used to describe 
wave functions of parent and daughter nuclei. Section Hill 
deals with the approximations employed and parameters 
used. In particular, we discuss the sensitivity of the cal- 
culated spectroscopic factor to the parameters defining 
the shifted Gaussian basis that is used to compute the 
normalization of the channel function. In Sec. IIVI we 
study the sensitivity of the reduced alpha width in ^^^Po 
on the choice of s.p. basis used. In Section |V] we discuss 
the absolute alpha-decay width of ^^^Po and in Sec. |Vl] 
we compare it with the absolute width of the superal- 
lowed alpha emitter ^""^Te. Finally, the main conclusions 
of this work are summarized in Sec. IVIII 



II. FORMALISM 

In this section, we discuss the R-matrix ([IJ and spec- 
troscopic factor ^ expressions for the decay width. 



The connection between the two formulations is given 
in Ref. Q. We also discuss the so-called delta- 
approximation for the formation amplitude. 



A. R-matrix expression for the decay width 

Within the R-matrix theory 0, the absolute 
width is given by Eq. ([iJ with Pl{R) being the barrier 
penetrability and jLiR) ~ the reduced width amplitude 
|27| . While both quantities strongly depend on the value 
of the channel radius R, the absolute width should be 
i?-independent. 

For Pl{R) we use the standard expression Q: 



Pl{R) 



kR 



\H+irj,kR)\^-' 



2m 



(3) 



ob- 



where k is given by the alpha energy Ea 
tained from the experimental Qa value by correcting for 

J^*^" is the reduced mass of 
alpha particle with 771^ being the mass of the daughter 

nucleus; H^{r],kR) is the outgoing spherical Coulomb- 

o 

is the Sommerfeld 



electron screening; fi 



Hankel function; and 77 = 
Coulomb parameter. 

The reduced width amplitude ^l{R) may be written 
in terms of the formation amplitude gL^R) 



y De w] 
I0: 



7i = 



9l{R), 



(4) 



with 



gdR) 



JM 



(5) 



where (j)a is the normalized wave function of the alpha 
particle with zero angular momentum, Ila/i, is the an- 
gular part of the center-of-mass (c.o.m.) motion of the 
alpha particle, ^'j^„^. is the wave function of the daugh- 
ter nucleus, and 'I'jm the wave function of the parent 
nucleus. The coordinates ^q. and are the intrinsic co- 
ordinates of the alpha particle and daughter nucleus, re- 
spectively. All wave functions are normalized in terms 
of the internal and c.o.m. coordinates [131 . By con- 
struction, the parent and daughter wave functions are 
antisymmetric. The antisymmetrization with respect to 
inter-fragment nucleons is done by means of the operator 

A. Its action can be approximated by means of a factor 

1/2 



("2" 



)( 



with N„ and Z„ 



being, re- 
spectively, the numbers of valence neutrons and protons 
in the parent nucleus. 

For the internal alpha-particle wave function we take 
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the standard Gaussian ansatz [1^, [20| : 

Xoo(cricr2) = [Xl/2(c^l)Xl/2('72)]oo, 



The parameter /3 = g^:5-=0.057fni ^ depends on the 

root-mean-square alpha radius r^^l.STfm (29| . 

The transformation between the intrinsic = 
{pi,P2,P3} and nucleonic {r^} (i=l,2,3,4) coordinates 
reads: 



Pi 

P2 
P3 



r-1 - r2 

r-3 - T-4 

(^"1 +^2) - (rs +r4) 



(7) 



and 



R 



ri+r2+ rs + 



(8) 



is the c.o.m. coordinate of alpha particle. Let us de- 
note the spherical components of intrinsic coordinates by 
Pi = {PiT^iiVi)- Assuming 6ji = ipn = 0, the nucleonic 
coordinates can be written as: 



4r^ 2 = 4i?^ 



pI 



^'3,4 



2p\ ± 2^2^3/31 cos 6^31 
+AR (^p3 COS 03 ± V2pi COS , 

4i?2 + p2 + 2pl T 2V2p3P2 cos 6^32 
-4i? (^P3 cos ^3 ± V2p2 COS 6*2) , 



(9) 



where 9ij — 6j — 9i, and 



cos Ui_2 



cos ^3^4 



2R + p3 COS 63 ± V^pi cos 6*1 

2R — p3 cos 63 ± V2p2 cos ^2 
2^4 ■ 



(10) 



This paper deals with g.s.— >g.s. alpha decays to the 
magic daughter nucleus. Assuming the seniority-zero 
wave function, the corresponding formation amplitude 
is [11,111 



FoiR) 



V8 
167r3/2 



where 



Iur„Up{R) ^ j dpidp2dp3(t){piP2P3) 



ri r2 



P,„(cos0i2) (12) 



^'ip (cos 6*34), 



with 9ij = 9j — 9i, V = {n,l,j}, and u^(r) being s.p. 
radial wave functions. The factor VS comes from the Ja- 
cobian of the transformation between the nucleonic co- 
ordinates {vi} and the internal and c.o.m. coordinates 
[ilM]. In Eq. (HI]) and in the following, the s.p. indices 
1,2 refer to neutrons while 3,4 refer to protons. The 
coefficients b^^_^^ are the shell- model four-particle wave 
function amplitudes. 



B. Delta-function approximation 

In the calculation of alpha-decay rates based on h.o. 
wave functions, it was noticed [2^ that the relative rates 
change little with the oscillator length 6h.o. of the basis. 
Using this argument, Mang proposed to take (3 ^ 1/6^ ^ . 
In this limit, the expression for the formation amplitude 
can be simplified (see also Ref. [sT']). In the literature, 
this is known as delta- function approximation 32]. 

In practice, one assumes that the alpha particle wave 
function is constant inside a small volume of radius Sq, = 
2.34 fm [32 j and zero outside. Within this approximation 
Pi = 0; hence, it immediatel y fo llows from Eqs. © that 
ri = r2 = r3 — r4 ^ R [S^, Isl], and the formation 
amplitude reduces to 



with 



(13) 



(14) 



where t — n,p. The correction factor depends on the 
relative angular momentum [32]: 



B^ = l- O.OUm^ + 1) 



C. Four-particle amplitudes 



(15) 



^ (-)'"+'^6,,.,,J„jp/,„,,^(ii'), (11) ^here 



For the g.s. alpha decay of ^^^Po and ^°'*Te, we are go- 
ing to assume that the four valence nucleons move around 
the rigid, doubly-magic core. The parent-nucleus wave 
function is approximated by a product of two-neutron 
and two-proton seniority-zero states: 

|*2n,0Q) ® |*2p,00), (16) 



l*J=0,A/=o) 



(17) 



|i^i/,00) = '"-^'"" lO^), and |0) = |0„) ® |0p) is the shell- 
model vacuum representing the ^^^Pb or ^""Sn g.s. wave 
function. The four-particle amplitudes b^„.Up in (|lip can 
thus be written in a separable form: 



(18) 
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D. Alpha decay spectroscopic factor 

Based on the general theoretical arguments [sl-flU , the 
absolute width can be expressed as a product of the 
alpha-particle spectroscopic factor and the single par- 
ticle width, see Eq. ([2]). The spectroscopic factor Sl 
contains information about the probability of forming an 
alpha cluster in the parent system. Since the alpha par- 
ticle, when formed, occupies the resonant state, the s.p. 
width can be obtained from the so-called current expres- 
sion ^,113: 

.h^ u'liR) ul{R) - u'l^iR) ul{R) 



where 



2^ 



J \ul{RW dR 



(19) 



where the Gamow function ui^{R) is obtained as a solu- 
tion of the Schrodinger equation with outgoing boundary 
condition. When the imaginary part of the complex en- 

fc2 7 2 

ergy eigenvalue £a = 2 small, which is always the 
case for the considered g.s. alpha emitters, one can ap- 
proximate (Uni) with [36|]: 



Y-sp 
^ L 



\ul{R)? 



\Ht{r^,kR)r 



(20) 



The s.p. width obtained in this way should be identical 
to the value — 2Im(fc) given by the imaginary part of the 
Gamow resonance energy, if the latter is computed with 
a sufficient precision. 

The conventional alpha spectroscopic factor as intro- 
duced in Ref. [8] is defined by 

Sl - |(^ [0o(ea) ^^fi^D) MR)],jM (21) 



where ipLAiiR) 



Ylm{R) represents the relative 



motion alpha particle with respect to the daughter, 
terms of the formation amplitude, Sl reads [sl. M |37| : 



In 



gl(R)R^dR. 



(22) 



Gl{R) = J U-^'\R,R') 9l{R')R'^ dR' (24) 

is the modified formation amplitude. The norm kernel 
Nl appearing in Eq. ([24]) is [43i] 



AfL{R,R') 
MRa~R) 



(25) 



D^ I A HR-a — R') 



i?'2 



[YL^f] 



3 \JI 



The presence of the norm kernel M effectively enhances 
the s pec troscopic factor by one-to-two orders of magni- 
tude [2l|, gO, Sl gl gi] . 

To compute M-^ {R,R'), we expand the eigenfunc- 
tions of the norm kernel in an orthonormalized shifted 
Gaussian basis (SGB) [H, 

FLiR,Rk) = y^(Np'^^) FL{R,Rk'), (26) 

— \ / kk' 



with Rk equidistant mesh points in the interval (0, -Rmax) 
and k — 1, . . . , M, where M is the dimension of the basis. 
The SGB is given by 



8/3' 



3/4 



FL{R,Rk) = 47r (^^j jL{-i^P' RRk), 

(27) 

while the SGB overlap {Np)kk' is given by 

{NF)kk' - J Fl{R,Rk)FLiR,Rk')R^dR (28) 
= 47re-'0'(^l+<'\^J^(-^4p'R,R,,). 

Using the SGB overlaps, the eigenvalue equation for the 
norm matrix can be expressed in the form: 



M 



(29) 



E. Modified spectroscopic factor 

Since the formation amplitude Eq. ([S]) represents the 
overlap of the parent wave function with the daughter- 
alpha product state, one would be tempted to associate 
it with the probability amplitude that in the parent wave 
function $jM alpha particle (/)q, and a daughter nu- 
cleus VP^L are at a distance R. The value oi Sl would 
then be associated with the total probability of formation 
of an alpha particle. However, the fundamental prob- 
lem with this interpretation is that the channel function 
■A r^(Ca) *f (C-d) V'l(-R)] jj^j is not properly normahzed 

|,[i3liil[3iSil. 

The properly defined spectroscopic factor (sometimes 
referred to as "the amount of clustering" ) j21|, 188,, ,42i - t45| 
is given by 



where 



Gl{R)R^dR, 



(23) 



— \ / kn 



L' [Np'^') „, (30) 

n k 



For /?' = 4/3, the core-projected norm Af^ in Eq. ((30)) 
reduces to a simple expression ^21, 43, 47]: 

^f^k> = ((^r'^^i^ir^^"))' ((^i^^'^i^i:^'^))' (31) 

where 

- J2 '5ii(4l^»y)(^ni,l4')(32) 
nljr^core 

with (^^(i?) = FL{R,RkW ^ 13) and i?„,,(i?) = u^i^/R 
are the radial s.p. wave functions of the core. 
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In terms of eigenstates of ([29]) . the spectral repre- 
sentation of the norm kernel can be written as: 

Af-'/'{R,R')= n-^^'u^*{R)^UR'), (33) 

(".^>"mi„) 

where the eigenfunctions (i?) of the norm kernel are 

M 

u^{R)^Y. ^^khiR^Rk), (34) 

k 

and rimin represents the usual cutoff on the eigenvalue 
of the norm kernel. The final expression for the modi- 
fied formation amplitude in the normalized SGB becomes 

Gl{R)= n-^'^ u';:{R) (35) 

(n..>n„i„) 

with 



TABLE II. The eigenstates (in MeV) of the s.p. Hamilto- 
nian of Table |I] for ^"^Pb calculated with the Gamow solver 
ANTI The positive-energy eigenvalues represent Gamow 
resonances; their imaginary energies reflect nonzero particle 
width. 



Orbit 


Neutrons 


Orbit 




Protons 




159/2 


-3.926 


Qhg/2 




-3.784 




Oill/2 


-2.797 


1/7/2 




-3.542 




2^5/2 


-2.072 


Oil3/2 




-1.844 




Oil5/2 


-1.883 


2p3/2 




-0.690 




3Si/2 


-1.438 


1/5/2 




-0.518 


11 


^"3/2 


n 781 

— u. * 01 


^Pl/2 


0.491 - 


- i0.200 X 10" 


I57/2 


-0.768 


I59/2 


4.028 


- i0.130 X 10" 


-7 


1/1-11/2 


2.251 -i0.026 


0«ll/2 


5.434 


- iO.992 X 10~ 


-8 


0jl3/2 


5.411 - i0.009 


0jl5/2 


5.960 


- j0.115 X 10" 


-7 






2d5/2 


6.748 


- i0.184 X 10" 


-2 






3Sl/2 


7.843 


- jO.367 X 10" 


-1 






157/2 


8.087 


- iO.898 X 10" 


-3 






2d3/2 


8.530 


- iO.284 X 10" 


-1 






lftll/2 


11.390 


- j0.215 X 10 


-1 






0il3/2 


15.086 


- iO.493 X 10 


-2 






l/l9/2 


15.964 - iO.393 





= / u^{R) gL{R) R^dR. (36) 



III. THE MODEL 
A. Single-particle space 

The s.p. space is spanned on resonant states of a 
WS-|-Coulomb average potential. The parameters of the 
s.p. Hamiltonian, namely the WS potential depth Vq, 
spin-orbit potential depth 14o, diffuseness a (= ago), ra- 
dius To (= ro^so), and the radius of the uniform charge 
distribution Tc defining the Coulomb potential are listed 
in Table ID The resulting neutron and proton s.p. ener- 



TABLE I. Parameters of the average WS Hamiltonian used in 
this work to compute s.p. neutron and proton states of ^"^^Pb 
and ^""Sn cores. 



Core 


T 


Vo 




a 










(MeV) 


(MeV) 


(fm) 


(fm) 


(fm) 


208pb 


n 


44.40 


16.5 


0.70 


1.27 




P 


66.04 


19.0 


0.75 


1.19 


1.27 


i««Sn 


n 
P 


51.60 
52.20 


11.3 
10.5 


0.70 
0.70 


1.27 
1.27 


1.27 



gies for 208pb and i°°Sn are given in Tables U and ITTIl 
respectively. The nucleus ^""'^Sb is proton- unbound; the 
values in Table IIIII are generally consistent with system- 
atics [48|. In particular, we predict a very small splitting 
between the 037/2 and 1^5/2 neutron shells outside ^™Sn, 
and a O57/2 g-s. in ^"^^Sn as suggested by recent experi- 
ment [il. 



TABLE III. Similar as in Table |lT] except for ^''''Sn. 



Orbit 


Neutrons 


Protons 


Off7/2 


-10.830 


2.669 -iO. 207 x 10" ' 


1^5/2 


-10.674 


2.869 - iO.963 x 10"^ 


2Si/2 


-9.074 


4.150 - j0595 X 10"^ 


1^3/2 


-8.927 


4.393 - i0.166 x 10"^ 


O/I11/2 


-5.793 


7.280 - iO.UO X 10"^ 


1/7/2 


-2.346 


9.649 - jO.452 


2p3/2 


-1.531 




2pi/2 


-0.912 




O/19/2 


-0.641 


12.012 - j0.0736 


1/5/2 


-0.171 




0«13/2 


3.254 - i0.132 x 10"^ 


15.572 -iiO.185 



B. Two-particle interaction 

The correlated two particle wave functions |5'2r.o) pT]) 
have been obtained using a separable two-body T = 1 
pairing interaction [sij : 

{lyu, 00\VW'i^', 00) = -Grfiiy, t)J{v', r), (37) 

where 

/(i.,r)-^-^(j,||ro|b.)/(^^,r). (38) 

In Eq. (l38l) we used the Condon-Shortley phase conven- 
tion for (ji/||lolbiy) and 

I{v,t) = / uUr)fr{r)dr. (39) 
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For the radial form factor fr (r) we took the derivative of 
the WS potential muhiphed by r: 



fr{r) 



(40) 



In the case of ^^^Po and ^""'Te the two-particle amph- 
tudes of Eq. (IT7)) were obtained exactly in the Tamm- 
Dancoff approximation (52l . [ssj : 



Xr. = Nn 



(41) 



where are s.p. energies, Eq is the correlated two- 
particle energy and A^o is the normalization constant 
fixed by the condition {XI f' — 1. 

The parameters R^r and a^i- defining the radial form 
factor (|40p for ^^°Pb and ^^°Po were chosen to repro- 
duce the wave functions used by Harada . Since such 
data are not available for ^''^Sn and ^°^Te, in this case 
we adopted the values of the WS potential for ^'"^Sn 
shown in Table HI The pairing strength Gr was adjusted 
to fit the experimental two-nucleon separation energies 
S2t through the dispersion relation 



2er. - et: 



(42) 



Since the proton- unbound nucleus ^''^Te is not known 
experimentally, for this system we adopted the value of 
S2p = — 2.14MeV obtained by extrapolating down from 



the heavier Te isotopes [5J|. This value is in reasonable 
agreement with recent phenomenological estimates p8| . 
Table IIVI lists the parameters of the residual interaction 
used in our study. 



TABLE IV. Parameters R^t and a^r of the residual interac- 
tion (|37() . The last column lists the value of S2t that has 
been used to constrain the pairing strength Gt for various 
configuration spaces considered. 



nucleus 


Rv (fm) 


ay (fm) 


52r (MeV) 




7.525 


0.70 


9.123 


210 


5.451 


0.75 


8.783 


lo^Sn 


5.895 


0.70 


24.3 


Te 


5.895 


0.70 


-2.14 



C. Configuration space 

To study the dependence of the formation amplitude 
on the size of valence space, and to compare with previous 
work, we considered several model spaces. Those used in 
the description of the alpha decay of ^^^Po are given in 
Table |Vl The model space MO contains only one valence 
shell. The space Ml contains one major shell, including 



the unusual-parity intruder orbit. The model space M2 
is that used by Harada [13]. The model space M3 is 
that of Glendenning and Harada [s^. Finally, M4 is the 
extended shell model space employed by Tonozuka and 
Arima. The model spaces used to describe ^^'^Te alpha 



TABLE V. Model spaces used in this work to describe ^^^Po 
alpha decay. 



Model 


Neutron States 


Proton States 


MO 


I39/2 


Ohg/2 


Ml 


lff9/2, Ojil/2, 2^5/2, Ojl5/2 
3Sl/2, 2(^3/2, 1^7/2 


Oha/2 , 1/7/2 , 0ji3/2 , 2p3/2 
1/5/2, 2pi/2 


M2 


lff9/2, Ojii/2, 2(^5/2 


OhQ/2, 1/7/2, 0ji3/2 


M3 


lff9/2, Oiii/2, 2^5/2, 0ji5/2 


O/19/2, 1/7/2, 0ii3/2 


M4 


lff9/2, Ojil/2, 2^5/2, 0ji5/2 
3Sl/2, 2d3/2, l(/7/2, 1^11/2 
0jl3/2 


Ohg/2, 1/7/2, 0ji3/2,2p3/2 
1/5/2, 2pi/2, lfl9/2,0iii/2 
0jl5/2,2d5/2,3Si/2, lff7/2 
2^3/2, l/lll/2, 0ji3/2, lhc,/2 



decay are shown in Table IVII Ml consists of one major 
shell, including the unusual-parity intruder orbit, while 
M4 consists of states with width less than 1 MeV. 



TABLE VI. Model spaces used in this work to describe " Te 
alpha decay. 



Model 


Neutron States 


Proton States 


Ml 


O57/2, ld5/2,2Si/2, 1(^3/2 

O/111/2 


0^7/2, 1^5/2, 2Si/2, 1^3/2 

Ohii/2 


M4 


O37/2, 1^5/2, 2Si/2, ld3/2 
0^11/2, 1/7/2, 2p3/2,2pi/2 
Oh,j/2, 1/5/2, Oil3/2 


0^7/2, 1^5/2, 2Si/2, 1^3/2 
Oftll/2, 1/7/2,0/19/2,0113/2 



D. Wave functions 

For the alpha formation amplitude in ^^^Po discussed 
in Sec. IIVI we considered the model spaces M2, M3 and 
M4. The wave function amplitudes in M2 were taken 
from Refs. [3, HB]- For calculations in M3, we took 
the T=l seniority-zero amplitudes of Ref. fssl] and renor- 
malized them accordingly. For calculations in the ex- 
tended space M4, we used the renormalized amplitudes 
of Ref. [l5|; here we retained only configurations hav- 
ing width smaller than 1 MeV. The comparison between 
^^^Po and ^°'*Te discussed in Sec. IVII was carried out in 
the model spaces Ml and M4. The corresponding wave 
functions were calculated in the two-particle approxima- 
tion described in Sec. IIIIBI except for ^^^Po in the M4 
model space, where Ref. ^53] was used instead. 



E. Penetration factor 

The s.p. alpha width Pq^ has been obtained from the 
current expression (I20p . The alpha-core potential was 
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assumed to be of a WS+Coulomb form with the param- 
eters of Ref. HI]: ro = Rc = 1.315 fm, a = 0.65 fm. The 
strength of the WS potential has been adjusted to repro- 
duce the measured Qa value corrected by the electron 
screening term [l^, l58i-,60.] : 



where 



1.4 



(eV). 



(43) 



(44) 



For 2i2po, Ea 8.785 MeV ^ and AE^^ ^ 31.8 kcV; 
hence, Qa = 8.986 MeV. The g.s. alpha decay of ^°^Te 
has not been observed. For that reason, we took the 
value (5^^=5.135 MeV extrapolated down from the bind- 
ing energy differences in ^°*Te (3.445 MeV) and ^°*^Te 
(4.290 MeV) [6lj. By adding the screening correction 
AEsc = 16.1 keV, we arrived at Q^ = 5.151 MeV. The 
resulting WS potential strength is Vq = 143.49 MeV for 
2i2po and 149.64 MeV for ^°^Te. 

The Gamow jvave functions were obtained by means of 
the code ANTI [50] . The complex energy of the metastable 
alpha state is = (8.986-i0.632 x IQ-^^) MeV for ^^"^Po 
and £a = (5.151 - i0.814 x 10"") MeV for ^^'^Te. The 
outgoing spherical Coulomb-Hankel function iJ+ was cal- 
culated using the code j63 |. 



F. Calculation of the spectroscopic factor 

The radial integration in the expressions for the spec- 
troscopic factor (1231) and the formation amplitude in the 
normalized SGB ((36)) have been carried out using 200 
Gauss-Legendre mesh-points with the maximum radius 
of 20 fm. 

The s.p. core wave functions entering Eq. (j32]) are 
those of the s.p. Hamiltonian of Table HI The radial 
mesh Rk defining the normalized SGB (|26l) was taken at 
equidistant points Rk = k AR. In order to determine the 
step AR we expanded the s.p. core states u{r) in the nor- 
malized SGB: u{r) = J^kLi '^fc [rFQ{r, Rk)]. Under the 
condition that Mdiff (^) = !"('") ~ ^'('") 1 1 < 0.005 fm"^/^ we 



found that 0.44 fm < AR < 0.57 fm and i?,„ax > 14 fm. 
For this range of AR and i?max the normalized SGB is 
orthonormal with an accuracy better than 10"^. To il- 
lustrate the quality of the resulting expansion, Fig. [T] 
shows UdiB{r) for the neutron core states in ^'^'^Pb. 

To calculate the modified formation amplitude G{R), 
one needs to determine the eigenvalue cutoff n^in- To 
this end, we show in Figs. [2] and |3] typical distribution of 
the eigenvalues of the norm kernel ((25)) for ^^^Po and 
^°*Te, respectively, for different values of AR. One may 
observe that a significant fraction of them accumulate at 
zero [53, 63]. To eliminate these spurious eigenvectors, we 
define the cutoff at the value where the eigenvalue distri- 
bution changes slope. For ^^^Po and -'^^'^Te this happens 
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FIG. 1. UdiB{r) — \u{r) — u{r)\ for the neutron core states 
in ^"^^Pb for AR = 0.5 fm, M = 30, and i?max = 15 fm. 
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FIG. 2. Eigenvalues of the norm kernel (|25|) for ^^^Po for 
fimax = 13 fm for different values of AR. 



at around 10^'^. Consequently, in our calculations, we 
adopt the cutoff value of n„iin — 0.001. 

The eigenfunctions (R) of the norm kernel ((34)) are 
orthonormal with an accuracy of lO"^'' for all eigenval- 
ues. The eigenfunctions with < nmin oscillate inside 
the nuclear volume and vanish outside the surface region. 
To further check the quality of (i?) we compute expres- 
sion (1351) by assuming nmin = and = I for all v. In 
this case, Eq. ((35)) formally reduces to g{R). Figure [4] 
shows g(R) for ^^^Po calculated in this way. The agree- 
ment with the original formation amplitude is excellent, 
except for a small deviation close to i? = and a small 
oscillation around and beyond the nuclear surface, which 
is not visible in the scale of Fig. Uj 

Next we study the sensitivity of S to the choice of 
i?max, AR, and nmin. For this analysis we relax the 
condition for Mdiff('') in order to access a wider range 
of AR. First, we study the sensitivity of iS as a function 
of i?max for various values of AR. Figure [5] shows the 
result for 0.53fm < AR < 0.59 fm for ^^^Po in the model 



8 



-01 
-03 

7 

^-05 
^-07 
-09 



-11 





1 ' 1 ' 1 ' 


1 1 1 1 1 1 1 1 1 1 1 1 1 






\\\> \ V - 




AR (fm) 


\v\\\ ^ 


r 


0.455 






0.476 


\aVVv i 




^ 0.500 






0.526 






^ 0.556 









' 



2 4 6 8 10 12 14 16 18 20 
eigenvalue number v 
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FIG. 5. Convergence of 5 for ^"Po (model space M4) as a 
function of i?max of the normalized SGB for different values 
of AR (with nmin = 0.001.) 
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FIG. 4. Formation amplitude g{R) for ^^^Po in the M4 
model space expanded in eigenfunctions of the norm kernel 
for AR = 0.500 fm and TJ^ax = 13 fm. 
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FIG. 6. Similar as in Fig. [5] but as a function of the step size 
AR for different values of i?max. 



G. Integral over intrinsic coordinates 



space M4 and nmin = 0.001. Except for a small value of 
AR = 0.53 fm, which does not produce stable results, a 
plateau in i?niax is reached around 14 fm. 

The dependence of 5 on AR displayed in Fig. [5] reflects 
the fact that for too small values of the step the basis 
functions become numerically linearly dependent, while 
for too large Ai?'s the basis cannot capture high Fourier 
components [H, H^lS^. Figure H] shows S for ^^^Po in 
the model space M4 and nmin = 0.001 as a function of 
AR. In general, appreciable oscillations of S can be seen 
except for the "safe" region 0.54 fm < AR < 0.59 fm, 
where results weakly depend on i?max- 

Finally, Fig. [7] shows the behavior of 5 as a function 
of the eigenvalue cutoff rimin for AR = 0.57 fm. The 

1/2 

cutoff used in Figs. [5] and [5] corresponds to n^-^^^^ = 
(0.001)-i/2 ~ 31.5. 



The multidimensional integral ()12p depends on the nu- 
cleonic coordinates, which are parametrized in terms of 
intrinsic variables through Eqs. (jH]) and (|10p . The inte- 
gration over (fii can easily be done analytically. Since the 
coordinates of particles 1 and 2 depend only on the rel- 
ative coordinates 1 and 3, and the particle coordinates 3 
and 4 depend only on the relative coordinates 2 and 3, 
one can greatly simplify the remaining six-dimensional 
integral by making first the integration over the relative 
coordinates 1 and 2 and then the integration over the 
coordinate 3: 



dps 



■ dp] 



■ dp-. 



(45) 



The integration has been carried out using the Gauss- 
Legendre quadrature using 10 points for the radial inte- 
grals and 8 points for the the angular coordinates. This 
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FIG. 7. Similar as in Fig. [5] but as a function of nmin for 
different values of -Rmax and AR — 0.57. 



guarantees the convergence up to the fourth significant 
digit. 



TABLE VII. Single-j alpha reduced width ratios at J? = 9.5 
fm. Shown are: rg of Ref. of Eq. (gS)), r of Eq. (|i7)) . 

and r of Ref. [65|1 . 



Orbital Jn 



0jl3/2 
Oill/2 
I99/2 
1/5/2 
2p3/2 



[32] 


rs 


r 


r [65] 


0.44 


0.46 


0.20 


0.10 


0.32 


0.34 


0.21 


0.08 


7.50 


7.50 


6.50 


3.73 


0.73 


0.74 


0.58 


0.55 


1.89 


1.89 


1.73 


1.89 



is seen that h.o. calculations underestimate WS values 
for high-j orbits by a factor two-to-three. 

It has been early recognized [1^ [s^ that the delta- 
function approximation overestimates the contributions 
of high-j orbitals. One can see it clearly by comparing the 
values of rg of Eq. with those of r (|Tf)) . i.e., the third 
and fourth columns of Table I VIII To cure this deficiency, 
a correction factor B^, was introduced [S^l in Eq. ([TJ) 
that depends on the relative angular momentum. 



IV. REDUCED WIDTH FOR 2i2po 
A. Single-j configuration 

Following Rasmussen [sS] , it is instructive to compute 
relative reduced widths assuming a pure single- shell 
model orbital assignment for the neutron pair, while the 
proton pair fills the Ohg /2 shell. For simplicity, the results 
are expressed relative to the ^^"Po reference (a neutron 
pair in 2^1/3). 

In the delta-function approximation of Sec. Ill B[ the 
ratio rs of the reduced widths is given by a simple ex- 
pression [32| : 



71 



rs 



2jn + 1 



'Pl/2 



(46) 



In a more general case expressed by Eq. (Ilip , the ratio r 
depends on the proton wave function: 



72pi/2>0/l9/2 



2j„ + 1 



hp^/2fihg/2{R) 



(47) 



Table IVIII compares the ratio rs given by Eq. us- 
ing the WS wave functions with that of Table I of Ras- 
mussen (32I based on the rounded square well potential 
of Blomqvist and Wahlborn [64] for several neutron con- 
figurations at i? = 9.5 fm. We find excellent agreement 
between these two calculations, and we checked that this 
agreement also holds for R — 9.0 fm. This is not surpris- 
ing as both calculations employ finite-depth potentials. 
The fourth column of Table IVlII displays the ratio r given 
by Eq. (|Tf| using the WS wave functions; they are com- 
pared with the h.o. values of Ref. [1^ (last column). It 



B. Eniiancement due to configuration mixing 

As was first shown by Harada the reduced width 
at the surface region is strongly enhanced by the config- 
uration mixing because contributions from various shell 
model orbits add coherently. To assess the effect of col- 
lective enhancement due to the configuration mixing, we 
carried out calculations in the M2 space. For i? = 8 fm, 
our WS calculations yield the enhancement factor of 
C = 8.5 with respect to the valence-shell configuration 
MO. This is to be compared with C = H obtained in the 
delta-function approximation; C = 10 obtained by Ras- 
mussen |32]; and C = 5.5 of Harada [H using h.o. wave 
functions. 

For the model space M3 of Glendenning and Harada 
[55| , obtained by adding the intruder neutron state 
to M2, we obtain C = 21. This should be compared 
with = 24 obtained in the delta-function approxima- 
tion and C = 30 obtained in Ref. [l^ (also within the 
delta-function approximation) using a fairly rich wave 
function that also includes proton-neutron correlations 
and J > two-particle couplings. It is worth noting that 
our enhancement is around 80% of that by Glendenning 
and Harada, and that the seniority-zero component in 
their wave function is also 80%. 



C. Extended shell model space 

Due to the strong collective enhancement of the re- 
duced width due to configuration mixing, it is important 
to consider extended shell-model space by taking into 
account higher- lying orbitals [l3|- For finite-depth shell- 
model potentials, such as the WS potential used in this 
study, this necessitates a proper treatment of the particle 
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continuum. An appropriate representation to deal with 
the continuum space is the complex Berggren ensemble 
representing bound and unbound s.p. states f66l.[67l[. 

Here we consider the large configuration space M4 of 
Tonozuka and Arima ^isj . i.e., all s.p. orbits up to TV = 7 
harmonic oscillator shell except for broad resonances with 
widths greater than 1 MeV. The shell-model amplitudes 
were taken from Ref. [l5| and renormalized to the re- 
duced model space. For the sake of comparison with 
Ref. 113, we consider the relative reduced width 



(48) 



where 7^ [R) = is the Wigner limit [68^ . 

Table IVIIII compares our WS results for 9"^ (R) with 
those of Ref. 15] obtained in the h.o. basis for several 
values of R. Generally, the reduced width obtained in 



TABLE VIII. The relative reduced width 9^ 
Ref. nSl and this work. 



obtained in 



Model space 


R (fm) 


Ref. [15] 


This work 


MO 


8.4 


6.3 X 10"'= 


0.60 X 10"'' 


M3 


8.5 


4.4 X 10~^ 


0.48 X 10"^ 


M4 


9.0 


2.9 X 10"* 


0.41 X 10"'' 



the WS model is about one order of magnitude smaller 
than that in the h.o. basis. This is because the h.o. 
basis knows nothing about the particle thresholds, and 
the radial behavior at large distances is solely determined 
by the oscillator length. For that reason, calculations 
based on the h.o. wave functions show large sensitivity 
to this parameter 
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FIG. 8. Formation amplitude g{R) for '^^'^Po obtained in this 
work in the model spaces MO to M4 as defined in table|Vl The 
imaginary part of the formation amplitude in M4 (dotted line) 
is also shown. 

The formation amplitude obtained in this work is 
shown in Fig. |S] for the configuration spaces MO, Ml, 



M2, M3, and M4. Compared with the formation ampli- 
tudes of Ref. [isj , the maximum of the formation ampli- 
tudes obtained in the WS model are significantly larger, 
and appear at lower values of R, than in the h.o. model. 
Also the overall shape of the formation amplitude is very 
different in the two cases. A characteristic two-humped 
shape of g{R) calculated in M4 resembles the formation 
amplitude G{R) obtained in Refs. H5i . i42i]. A similar re- 
sult was also obtained in Refs. (20ll2l|. It is indeed inter- 
esting to see that a two-humped behavior of the forma- 
tion amplitude for ^^^Po has been obtained by consider- 
ing large configuration space and the Berggren ensemble 
of the WS potential. 

Figure |8] also shows that the formation amplitude in 
the M4 model space has a small imaginary part. This 
is because our calculations are carried out in the pole 
approximation that ignores the non-resonant continuum 
[67l [70l - l72| . This spurious component of g{R) results in a 
very small imaginary contribution to the reduced width, 
which can be safely neglected considering the expected 
accuracy of our model. 



V. ABSOLUTE ALPHA-DECAY WIDTH OF 

212p^ 



The g.s. alpha-decay width of ^^^Po has been deter- 
mined in the seniority-zero approximation using three 
different models spaces listed in Table El MO, M3, and 
M4. The corresponding four-particle shell-model wave 
function contains one configuration in the MO space, 12 
configurations in M3, and 144 seniority zero configura- 
tions in M4. 

The absolute width from Eq. ([T|) should not depend 
on the channel radius R. However, in R-matrix stud- 
ies involving approximations, such as the one-channel 
R-matrix treatment, this condition cannot be met p^ . 
Therefore, in practical calculations, in which the depen- 
dence of Tl on R around the nuclear surface is small rel- 
ative to the appreciable i?-dependence of the formation 
amplitude, one is trying to meet the plateau condition for 
r L (R) in which the absolute width varies weakly around 
the nuclear surface (T^I . Figure [HI shows the dependence 
of the R-matrix width ([T]) on the channel radius. It is 
seen that the plateau condition is met only in the case of 
the extended configuration space M4 involving particle 
continuum. Here, we find a fairly weak variation of F(i?) 
between 7fm and 11 fm. 

As seen in Fig. [9l and discussed in Sec. IIVBI and 
Refs. [15, 69, 74], the width strongly increases with the 
size of the shell-model space. Indeed, in the surface re- 
gion, F(i?) obtained in M3 shows an enhancement ~15 
with respect to MO, and in the extended space M4 the 
enhancement is ^^260. Compared to experimental value, 
however, the width obtained in M4 is still 600 times 
smaller than the experimental value Fcxp — 0.153 x 10"^'' 
MeV H. 

A further enhancement in the reduced width is due to 
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FIG. 9. Dependence of the absolute alpha-decay width ((T]) 
of ■^^■^Po on the R-matrix channel radius R for three different 
model spaces MO, M3 and M4. 



the antisymmetrization and normalization of the channel 
decay [H, . This is achieved by replacing the stan- 
dard formation amplitude g{R) with the modified forma- 
tion amplitude G{R) of Eq. i^. Figure [13 shows G{R) 
calculated in the M4 model space with AR — 0.56 fm, 
-Rmax = 11-76 fm (M = 21) and n^^n = 0.001. A smah 
oscillation at the tail of G{R) can be seen. The amplitude 
of this oscillation, around the asymptotic behavior given 
by Hq^t], kR), varies very little with_Rmax for this value 
of A^. As discussed in, e.g., 0, IH, S^l, the behavior of 
g{R) and G{R) is generally very different. This can be 
seen by comparing Figs. [51 and [TUl 
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FIG. 10. Modified formation amplitude G{R) of Eq. ([Hi in 
the extended model space M4 with rimin = 0.001, AR = 0.56 
fm and JJmax — 11.76. Unlike g{R), G{R) properly accounts 
for the normalization and antisymmetrization of the decay 
channel. The asymptotic behavior of G{R) is given by the 
Coulomb-Hankel function at the alpha-decay energy Qa ~ 
8.986 MeV (dashed line). 

The absolute alpha-decay width obtained by using the 
R-matrix expression ([1} with the formation amplitude 



G{R) of Fig.[TO]is shown in Fig.[TT] There appears a small 
plateau in the region of nuclear surface that corresponds 
to F « 0.0042 X 10^1^ MeV. This value is - 36 times 
smaller than Fcxp- At larger distances R > 9fm, the 
result is affected by spurious oscillations of G{R) around 
Hq {ri,kR), i.e., it is quite unreliable. 
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FIG. 11. Absolute width from R-matrix expression ([T]) cal- 
culated in the M4 model space using the modified formation 
amplitude G{R) of Fig. 1101 At 7? > 9fm, the result obtained 
by assuming G{R) oc Hg{T}, kR) is marked by a dotted line. 

The absolute width can also be obtained from expres- 
sion which involves the alpha-particle spectroscopic 
factor S and the s.p. decay width. Figure [1^ shows the 
result of the current expression (I^Ul) for F^^ as a function 
of the channel radius. As discussed in Ref. [36[, F'*^ cal- 
culated this way should be independent of i? if ii is large 
enough. This is precisely what is seen in Fig.[T2 the s.p. 
width converges beyond the range of the WS potential to 
VP = 0.1247 X 10-12 ]y[gv, which is indeed very close to 
the value of -2Im{£a) = 0.1265 x lO^^^ y^^y given by 
the imaginary part of the Gamow resonance. 

Using the modified formation amplitude G{R) of 
Fig. [nl we compute the spectroscopic factor S — 0.011, 
which - combined with the value of F''^ above - yields 
F = 0.14 X 10"" MeV. Using Ai? = 0.55 fm we obtain 
S = 0.0080 and F = 0.10 x 10"" MeV. Both these values 
are close to Fexp = 0.153 x lO^^^MeV. 



VI. COMPARISON BETWEEN 
GROUND-STATE ALPHA DECAY OF ^^^Po AND 

104 

To compare absolute widths of ^-'^^Po and -'^^^Te in a 
consistent way, we consider similar Ml and M4 model 
spaces for both nuclei. The norm kernel eigenvalues 
do not depend on the model space in which g{R) is cal- 
culated, so we take the cutoff n„iin = 0.001. 

Let us begin with ^^^Po by making a convergence anal- 
ysis of S in the Ml model space as a function of Ai? and 
i?max (as in Fig. O. For Ai? = 0.53, 0.54, 0.55, and 
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FIG. 12. Single particle width of '^^^Po from current expres- 
sion (f20)l . 

0.56 fm, we found S = 0.0041, 0.0011, 0.00030, 0.00032, 
respectively. The resulting converged value S — 0.0003 
is too small, as expected from Fig. El This deficiency 
is related to the poor quality of the interaction used to 
describe ^^^Po in Ml. To better understand this fact, let 
us take a look of the spectroscopic factor in terms of the 
spectral representation of the norm kernel, 

5 = (49) 

where the sum is truncated by the condition n^, > rimin- 
The summation range and eigenvalues Ui, are the same for 
Ml and M4; the only difference comes from g^. Because 
of the rapid oscillation of the eigenfunctions inside the 
nucleus, only the eigenfunctions which are peaked at and 
beyond the nuclear surface will contribute significantly 
to the sum. But - because g{R) in Ml is small in the 
surface region - the overlap with those eigenfunctions is 
small, and this gives rise to a very reduced value of S. 

By making a similar analysis for ^°^Te in Ml, we 
found S = 0.067, 0.024, 0.0066, and 0.00046 for AR = 
0.53, 0.54, 0.55, and 0.56 fm, respectively. In the model 
space M4 we found S = 0.21, 0.088, 0.032, and 0.0051 
for the same values of AR. Clearly, the convergence in 
S has not been achieved for ^°*Te. We would like to at- 
tribute this to the impact of the proton continuum on 
g^, which results in increased oscillations of G{R) in the 
surface area. Table IIXI compares the values of S and 
the corresponding absolute widths for ^^^Po and ^°'*Te 
at AR = 0.56 fm. (The single particle width for Te is 
T'P = 0.162 X 10-12 MeV.) 

It is interesting to compare our current results for ^''^Te 
with the estimates of phenomenological alpha-decay 
models based on semi-classical approximation 75 77). 
The assumed large value of Qa — 6.12 MeV in Rcf. [75| 
results in a very short half-life of 7 x 10-^^ sec. The alpha- 
decay energies of 5.05 MeV ^ and 5.42±0.07MeV ^ 
result in T1/2 ~ 10"'' sec and ^ 5 x 10~^sec, respec- 
tively, and these estimates are not inconsistent with our 



TABLE IX. Alpha decay spectroscopic factor and absolute 
width for ^^^Po and ^"^Te computed in the configuration 
spaces Ml and M4, with rimin = 0.001, AR = 0.56 and 
i?max = 11.76 fm(M = 21). 



Model Space 




5 


r x 10" 


MeV 




212p^ 


104^g 


212p^ 




Ml 


0.00032 


0.00046 


0.0040 


0.0075 


M4 


0.011 


0.0051 


0.14 


0.083 



value (M4 model space) T1/2 = 5.5 x 10"'' sec {Qa = 
5.151 MeV). As the value of Qa in ^'^''Te is very uncer- 
tain, we show in Fig. [13] the absolute width and half- life 
Ti/2 as a function of Qa for the model space M4. 
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4.9 5 5.1 5.2 5.3 
Q„(MeV) 

FIG. 13. Ground-state alpha-decay width (left scale) and 
half-life (right scale) in ^''''Te as functions of the decay energy. 

Our predicted spectroscopic factors in M4 for ^"^Te 
and -^i-^Po are about 0.5% and 1%, respectively. As men- 
tioned above, a fairly small value of S in ^'''^Te could be 
a consequence of the proximity of the proton continuum. 
Indeed, all the valence proton shells are resonances. The 
small value of S in ^^"'Te could also be attributed to the 
poor quality of the valence interaction assumed, and the 
neglect of the T = force. The effect of the proton- 
neutron interaction was examined in, e.g., Refs. |l4l. pTSj 
for -^i-^Po and was found to be minor due to the fact that 
neutrons and protons in ^^^Po occupy different shells. 
This is no longer true in the N = Z nucleus ^^^Te, in 
which the major enhancement of S is expected due to 
T = correlations. Therefore, our predictions for S and 
r in i"^Te given in Table IIXI should be considered as a 
very conservative lower limit. 

VII. CONCLUSIONS 

The g.s. alpha decay of ^^^Po has been studied 
within the complex-energy shell model framework with 
the Berggren ensemble of the average Woods-Saxon po- 
tential. We applied the pole approximation by consider- 
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ing s.p. resonant states only. The overlap integral involv- 
ing alpha-cluster nucleons was computed exactly, without 
resorting to the delta-function approximation. We con- 
sidered the large valence space of Tonozuka and Arima 
that is necessary to produce the collective enhancement 
of the formation amplitude. 

The absolute alpha-decay width was computed using 
the reduced width obtained in the framework of the R- 
matrix theory and also from the alpha spectroscopic fac- 
tor. The latter approach yielded results consistent with 
experimental value, but only after considering the anti- 
symmctrization and normalization of the decay channel 
wave function. The R-matrix estimate underestimates 
the experimental width by a factor of ^ 36. The R- 
matrix expression depends on the asymptotic value of the 
formation amplitude that is very sensitive to the size of 
the configuration space. On the other hand, the reaction- 
theory expression ^ involves the spectroscopic factor ~ 
an integral quantity that depends less on the size of the 
basis used. It is very encouraging to see that a reasonable 
agreement with the experimental width of ^^^Po has been 
obtained without explicitly considering the alpha-cluster 
component in the wave function of the parent nucleus. 
In this context, we believe that the improved treatment 
of the particle continuum has been essential. 

We have also provided an estimate of the alpha-decay 
rate in ^°^Te. Unfortunately, due to the fact that the 
valence proton shells in this nucleus lie in the continuum, 
no fully convergent result has been achieved. We hope 
to improve the situation in the future by inclusion of 
the non-resonant continuum space that will remove some 
of the undesired oscillations in G{R) at large distances. 
In addition, since the residual interaction employed in 
our work neglects the proton-neutron components, and 



the wave function has a seniority-zero character based 
on T = 1 nucleonic pairs, the predicted alpha width in 
this N = Z nucleus should be viewed as a conservative 
low limit. Indeed, the inclusion of T = correlations is 
expected to increase the value of F significantly. 

The calculations presented in this study should be 
considered as an important step towards an improved 
microscopic understanding of the alpha-decay process. 
Still, as this work demonstrates, further improvements 
are needed. The neglect of the non-resonant continuum, 
i.e., complex-energy scattering states in the Berggren en- 
semble, slightly violates the completeness relation at a 
one-body level. This results in small imaginary contribu- 
tions to spectroscopic factors and reduced widths, and - 
most importantly - can affect the behavior of formation 
amplitudes at very large distances. The second crucial 
development will be the use of large-scale shell model 
calculations, including realistic T = and T = 1 in- 
teractions, to compute wave function amplitudes. This 
will enable us to provide a more meaningful estimate of 
^°^Te alpha decay rate. The work in both directions is 
underway. 
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